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Thin nanomagnetic rings have generated interest for fundamental studies of magnetization reversal 
and also for their potential in various applications, particularly as magnetic memories. They are 
a rare example of a geometry in which an analytical solution for the rate of thermally induced 
magnetic reversal has been determined, in an approximation whose errors can be estimated and 
bounded. In this work, numerical simulations of soft ferromagnetic rings are used to explore aspects 
of the analytical solution. The evolution of the energy near the transition states confirms that, 
consistent with analytical predictions, thermally induced magnetization reversal can have one of 
two intermediate states: either constant or soliton-like saddle configurations, depending on ring size 
and externally applied magnetic field. The results confirm analytical predictions of a transition in 
thermally activated reversal behavior as magnetic field is varied at constant ring size. Simulations 
also show that the analytic one dimensional model continues to hold even for wide rings. 



I. INTRODUCTION 

The study of magnetization reversal in meso- and 
nanoscale ferromagnets is motivated by its importance in 
information storage applications. Ring-shaped thin film 
nanomagnets are particularly interesting because the ab- 
sence of sharp edges inhibits nucleation which may pre- 
cipitate undesired reversals. In addition, the geometry 
allows an analytic solution - rare in micromagnetics - for 
thermally induced reversal [ij . 

The lowest energy configurations of a soft ferromag- 
netic ring are a clockwise and counterclockwise circula- 
tion of the magnetization. The field generated by a cur- 
rent / flowing through the axis of the ring (z-direction) 
will lift the degeneracy between the two orientations. 
For an effectively one-dimensional ring (ring width small 
compared to the average radius, so that external mag- 
netic field strength within the ring is approximately con- 
stant), the reversal from the metastable state to the sta- 
ble state can occur through either of two processes: a 
global rotation of the magnetization (which we will call 
the constant saddle state), or reversal initiated within 
a localized region and subsequent expansion (instanton 
saddle) Q. Which of these processes occurs depends on 
both ring size and applied magnetic field. Using micro- 
magnetic simulations in OOMMF 0] we have confirmed 
this general picture, even for more general two dimen- 
sional rings. 

The model introduced in [3| predicted that the mag- 





Figure 1: Magnetization configurations for h = 0.2. a) Con- 
stant saddle, b) Instanton Saddle at I = 12, c) Instanton 
saddle at Z = 60. 



netization will be constrained to the film plane (m^ = 0) 
and the transition states will depend on both ring size 
and magnetic field. Important field and length scales in 
the problem are: 



(1) 
(2) 



Here Ms is the saturation magnetization, t is the ring 
thickness, Ai? is the ring width, R is the average radius, A 
is the exchange length. Hp is the external magnetic field, 
and He is the field at which the metastable configuration 
becomes unstable. Fori < 27rVl — h"^ the theory predicts 
a constant saddle, as shown in Fig. la, whereas for / > 
27rVl — /i^, it predicts an instanton saddle (Figs.[T|3, c). 
Both of these saddle configurations are described by a 
fimction (ph,i{d) [if. 



II. MICROMAGNETIC SIMULATIONS 

We simulated the magnetization dynamics of rings of 
mean radius R = 200 nm and thickness t ~ 2 nm, for ma- 
terials of exchange length A ranging from 4 to 40 nm and 
constant saturation magnetization Mg = 8x 10^ A/m (as 
in permalloy). The applied field has radial dependence 

^('') ~ 2TTr^ ~ where He is the magnitude of the 

magnetic field generated by the current / at the mean 
radius of the annulus. We studied three different regimes 
for I: (1) I < 2tt, where the transition state is predicted 
in [l| to be the constant saddle configuration (Fig. la); 
(2) I greater than but close to 27r and (3) Z >> 27r. The 
latter regime is relevant to recent experiments on nano- 
magnetic annuli [1, i, d, H, 0, Hi ■ 

The magnitude of the exchange stiffness (and there- 
fore exchange length) was varied so that I lay in all three 
regimes: respectively, I ~ 6, 12, and 60. Two different 
ring widths were studied: Ai? = 40 nm {He — 73.9 mT) 
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and AR = 100 nm {He = 29.5 mT). For different values 
of the parameter h the initial configurations (la), (lb) 
and (Ic) of the magnetization are allowed to evolve fol- 
lowing the Landau-Lifshitz-Gilbert equations: 



h,=0.21 
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where 7 is the gyromagnetic ratio and i?efF is the effective 
magnetic field (i.e., including both the external field He 
and internal fields generated by the ring magnetization) . 
The damping coefficient a can be set to one without loss 
of generality. 

We first conducted tests to determine how well the 
analytical solution presented in [l[ approximated the nu- 
merically determined two-dimensional saddle state in a 
narrow ring. We did this by initializing the spin con- 
figuration at time i = using the instanton configura- 
tion given by Eqs. (8-11) of |l|] with a fixed parameter 
h. We then allowed the system to evolve in time ac- 
cording to ([3]) with different values of the external field 
(of course, the external field is fixed during any particu- 
lar run). A saddle state for a particular external field ht 
must have the property that, for dynamics determined by 
slightly higher fields, the system falls to the stable state; 
for slightly lower fields, to the mctastable state (or vice- 
versa). Therefore, for a fixed initial spin configuration 
defined by h we varied the field determining the dynam- 
ical evolution until a "transition field" ht was found: for 
external fields less than ht the system would evolve to the 
metastable configuration while for those greater than ht 
it would evolve to the stable configuration. In this way 
the saddle configuration can be numerically determined. 
Values of ht were obtained with a numerical uncertainty 
oi5ht = 6 X 10-3. 



III. RESULTS AND DISCUSSION 

The above discussion implies that the analytical solu- 
tion is a good approximation to the actual configuration 
if the field h that determines the initial instanton con- 
figuration is close to the numerically dtcrmined ht- The 
results are shown in Fig. 2, which shows the the time 
evolution of the energy. Typically, after a short tran- 
sient the system arrives at a configuration for which the 
energy stays almost constant in time before decaying to 
either of the stable states. As h gets closer to ht, the 
system maintains a constant energy for longer times. We 
interpret this behavior as evidence that the analytically 
determined instanton predicted in [11] is indeed a close 
approximation to the actual saddle configuration. It can 
also be observed from Fig. 2 that at ht = 0.21 the energy 
of the constant saddle is higher than that of the instanton 
saddle, as predicted in [l[. The slope of the curve during 
the subsequent relaxation depends on how fast the rever- 
sal propagates along the ring: as A decreases the reversal 
time increases because both the domain wall width de- 
creases and the effective circumference I increases. This 
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Figure 2: Time evolution of the total energy in a field ht = 
0.21, for the instanton {h = {0.225, 0.25}) and constant saddle 
configurations (h = {0.3, 0.325}) when 1=12 and Ai?=40 nm. 
These transition states are chosen to bracket the saddle state 
at the field of the simulation - the numerically determined 
ht is not the same as the parameter h used to generate the 
saddle configurations from Eqs. (8-11) of [l[ (see text for a 
fuller explanation). The arrows represent the magnitudes of 
the activation energy AE{ht = 0.21) for each of the saddle 
states. 



reduces the magnitude of the slope of the curve E{t) dur- 
ing the reversal. 

The reversal activation energy AE is summarized in 
Fig. [3] as a function of AR and A. Every datum point 
corresponds to a step of Ah = 0.1 in the parameter of 
each of the saddle configurations. For fixed I, the ener- 
gies of the instanton and saddle configurations cross at 



he = y 1 — [^) ;at higher fields the constant saddle is 
the lowest energy transition state. As predicted in [l|, 
the instanton saddle has lower energy than the constant 
saddle for h < he, while for larger fields the constant sad- 
dle has a lower energy (Fig. 3a, h = 0.57). One should 
observe, however, that the predicted external switching 
field is larger than the numerically one. This is probably 
due to round-off from the simulations, but more impor- 
tantly from the fact that the theoretical result is for a 
one-dimensional geometry whereas our simulations run 
in a (more realistic) two-dimensional mesh. So it's sur- 
prising that the agreement between the theoretical and 
numerically observed values for the switching field im- 
proves when the width AR is increased from 40 nm to 
100 nm; cf. Figs. 3a-3c and 3b-3d. We believe this to be 
a consequence of the increased number of active cells in 
the simulation for the wider ring {AR — 100 nm); it im- 
plies that the Id solution is a good approximation even 
for fairly wide two-dimensional ring geometries. 

The theory presented in [l| predicts that the difference 
between the activation energy of the constant saddle and 
the instanton saddle increases as A decreases. This is 
also confirmed by our numerical simulations (compare 
Fig. 3a to 3b and 3c to 3d). This decrease is due to 
fact that the activation energy for the constant saddle 
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Figure 3: Energy barrier dependence on ht for both saddle 
configurations in different regimes. Activation energies for 
I = 12, AR = 40 nm at ht = 0.21 are sfiown for comparison 
witli Fig. 1. Data points represent h = {0.1, . . . , 0.9}. For 
comparison, theoretical lines are plotted for the activation 
energy. The crossing of the analytical curves in (b) and (d) 
is not visible due to the graphical resolution. 

is independent of A while the activation energy for the 
instanton (whose size is order A) is roughly proportional 
to A. This latter dependence of activation energy on A is 
in turn understood by observing that the total energy of 
the ring is the result of three contributions: Zeeman, ex- 
change and magnetostatic energy. The constant saddle, 
metastable and stable states all have the same exchange 
energy, so the exchange energy does not contribute to 
the energy barrier of the constant saddle. Both remain- 
ing terms (magnetostatic and Zeeman energy) are inde- 
pendent of A, so the constant saddle activation energy 
does not change with A. On the other hand, in most of 



the ring, the metastable and instanton configurations are 
almost identical, so that their energy difference is the re- 
sult of a fluctuation in a small region of the ring, whose 
size is of order A. The barrier (both theoretical and ex- 
perimental) for the instanton saddle is therefore smaller 
in Figs. 3b and 3d than in Figs. 3a and 3c. 



IV. CONCLUSION 

In these studies we have conducted numerical micro- 
magnetic simulations of the transition states for mag- 
netization reversal in thin annular nanomagnets to test 
the theory presented in [ij. Numerical calculations are 
consistent with the prediction that in the parameter 
space defined by {h, I) there is a region for which a non- 
constant, soliton-likc transition state (the "instanton" 
state) exists and has lower energy than the constant tran- 
sition state. An important assumption made in [l| is that 
the magnetostatic contribution is dominated by the sur- 
face magnetostatic effect. This is supported by the above 
simulations. One would expect that by increasing the 
ring width the bulk magnetostatic energy would increase 
in magnitude and eventually become nonnegiigible, so 
that other transitions paths become more probable (e.g., 
vortex-like states). This should eventually occur; how- 
ever, up to the appreciable AR/R values studied here, 
the approximations made in [l| remain valid. 
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